Orbitally relieved magnetic frustration in NaVC>2 
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The magnetic properties of NaV02 are investigated using full-potential linearized aug- 
mented plane wave method. We perform calculations for three structures. For the rhombo- 
hedral structure at 100 K, the t2 9 orbitals of V ions are split into upper a± g and lower e' g 
orbitals by a trigonal distortion of compression. For the monoclinic structure at 91.5 K, the 
system behaves like a frustrated spin lattice with spatially anisotropic exchange interactions. 
For another monoclinic structure at 20 K, the magnetic frustration is relieved by a lattice dis- 
tortion which is driven by a certain orbital ordering, and the long-range magnetic ordering 
is thus formed. Moreover, the small magnetic moment originates from the compensation of 
orbital moment for the spin moment. 
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I. INTRODUCTION 



The compounds with a common chemical formula AT0 2 (A=Na or Li, T=3d transition metals) 
have been attracting a lot of attention for their large variety and richness in physical phenomenajl , 
QflflBQ. Indeed, the discovery of supercondctivity in the Na . 35 Co0 2 -l 3H 2 O|10 and the appli- 
cation of LiCo0 2 in rechargeable Li batteries have accelerated investigations on their fundamental 
physics. In addition, NaMn0 2 undergoes a structural phase transition at 45 K to a long-range 
ordered antiferromagnetic (AFM) ground state[2], while NaNi0 2 exhibits ferromagnetic (FM) 
coupling in Ni-Ni plane below transition temperature yfl. Furthermore, a well-known member of 
this group, LiNi0 2 , has no long-range magnetic ordering even at low temperaturesjULM, III |9J]. 
The orbital frustration has been used to explain the absence of magnetic ordering |4|, [I ISfl, and a 
local ordering of Ni 3+ Jahn-Teller (J-T) orbitals is also proposed to be responsible for the com- 
plex magnetic properties [9]. Therefore, its controversial magnetic properties have been attracting 
considerable interest. 

These various phenomena always relate to the quasi-two-dimensional (2D) triangular lattice 
formed by T cations. Such a triangular lattice may lead to magnetic frustration, since all nearest- 

Finn, 

neighbor AFM interactions can not be satisfied simultaneously [2, 5, 61 . Nevertheless, the magnetic 



frustration can be relieved by a certain orbital ordering (QQ) [|1(a IllL 11211 . such as in LiV0 2 ll3j] 



and NaV0 2 dl4|]. LiV0 2 has been found to form a spin- singlet phase with corresponding 00 at 
low temperatures J5, 13]. Whereas, its sister compound NaV0 2 displays very different behaviours. 



Recently, Onoda et al. ||15|1 have revealed a superparamagnetic state driven by the short-range 
ordered spin-1 (the total spin in one trimer S=l) trimerization[16] below the transition temperature 
(r=98 K). However, McQueen et al. have reported two successive 00 transitions in NaV0 2 [14]. 
At 98 K, the system undergoes a continuous phase transition from a rhombohedral (R-3m) phase 
to a monoclinic (C2/m) one, corresponding to the proposed 00 of one electron per V 3+ . Below 
93 K, the system undergoes a discontinuous phase transition to another monoclinic (C2/m) phase, 
consistent with the proposed 00 of two electron per V 3+ . In addition, a long-range ordered AFM 
state is formed at low temperatures, while the magnetic moment observed in the ordered phase is 
about 0.98 \lb, much smaller than the expected value (2 /i#). The controversial magnetic states 
below 98 K obtained by these two groups bring us interests, and the puzzling magnetic moment 
deserves to be explored. There are no theoretical reports yet to the best of our knowledge, therefore 
we expect to understand its magnetic properties at low temperatures upon our theoretical efforts. 
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In the present work, we have performed first-principles calculations to investigate the electronic 
structures of NaV0 2 , further to reveal the most possible orbital and magnetic ordering, and to 
explore the origin of small magnetic moment observed at low temperatures. Partially in agreement 
with the experimental findings ilM, the 00, accompanied by a long-range magnetic ordering, is 
found for the second monoclinic structure. And the observed magnetic moment can be explained 
by including the spin-orbit coupling (SOC) interactions. 

This paper is organized as follows. The crystal structure and computational details are de- 
scribed briefly in Sec. II. And the results and discussions are presented in Sec. III. Finally, a brief 
conclusion is summarized in Sec. IV. 



II. CRYSTAL STRUCTURE AND COMPUTATIONAL DETAILS 



The lattice parameters provided by McQueen et al. are listed in Table I. The structure of 
NaV02 is composed of 2D triangular-lattice VO2 layers of edge sharing V06 octahedra separated 
by sodium ions, which is rhombohedral (R-3m) at relatively high temperature (HT) (T>9& K) and 
monoclinic (C2/m) in both intermediate temperature (IT) (91.5 K<T<98 K) and low temperature 
(LT) (7/<91.5 K) phases[14]. At 100 K, the V-0 distances are nearly 2.04 A, but the O-V-0 angle 
a is only 85.53° (Fig. 1). With lowering the temperature, a further reduces to 85.42° at 91.5 K 
and 85.13° at 20 K. Therefore, even in the HT phase, V06 octahedra in NaV02 have been different 
from the regular ones under a trigonal distortion of compression along the threefold (111) axis[14], 
which induces the lowering of the Oh local symmetry to D 3d . In V0 2 layers, the V-V geometry is 
built by two long (3.00781 A) and four short (2.99316 A) bonds in the IT phase, and reversely is 
by two short (2.97551 A) and four long (3.00526 A) ones in the LT phase llnll . The interlayer V-V 
distance is about 5.6 A. 

All the calculations were performed by using the standard full-potential linearized augmented 
plane wave code WIEN2k [|18|l . The muffin-tin sphere radii of 2.22, 2.00 and 1.77 a.u. were cho- 
sen for the Na, V and O atoms, respectively. The cutoff parameter R rnt K max was chosen to be 7.0 
and 100 ^-points were used over the first Brillouin zone. The local spin density approximation 
(LSDA) of Perdew and Wang lll9n was used for the exchange and correlation potential. In order to 
take the strong-correlated nature of 3d electrons into account explicitly, we performed LSDA+U 
calculations [20], where U e jf=U-J (U and J are on-site Coulomb and exchange interaction respec- 
tively) was used instead of C/Q2JJ] . And the orbital-dependent potential has the form of AVfll 



= -U e ff(n a - [22b, where h a is the orbital occupation matrix of sp 



spin a. This type of double- 



2411 . For NaV0 2 , we used 



counting correction (DCC) has been called the fully localized limit [|23L 
U e ff=3.6 eV which has been used in its sister compound LiV0 2 |15j]. Note also that the conclusion 
made in this paper is not affected for U e ff =2-6 eV|25|]. To explore the origin of small magnetic 



moment observed in the LT phase, we performed LSDA+SOC+U calculations, where the SOC is 



included by the second- variational method with scalar relativistic wave functions H18H . The easy 
magnetization direction was set along (1 10) direction (short V-V bonds in the LT phase) observed 
in the experiment ll 1411 . 

In order to investigate different magnetic patterns, 2x2x2 supercell was used in our calcula- 
tions. We took into account two AFM structures in V-V plane as described in Fig. 2(a): I- type 
antiferromagnetism (Fig. 2(a)(i)) is AFM exchange along the (010) and (110) directions with FM 
exchange along the (100) direction, II-type antiferromagnetism (Fig. 2(a)(ii)) is AFM exchange 
along the (100) and (010) directions with FM exchange along the (110) direction. Totally there 
were five possible magnetic configurations in our calculations for the IT and LT phases (Fig. 2(b)): 
ferromagnetic (FM), C-AFI (I-type antiferromagnetism in plane, FM stacking), C-AFII (II-type 
antiferromagnetism in plane, FM stacking), G-AFI (I-type antiferromagnetism in plane, AFM 
stacking), G-AFII (II-type antiferromagnetism in plane, AFM stacking). 

III. RESULTS AND DISCUSSIONS 



A. HT phase 

As the paramagnetic behavior of NaV0 2 has been determined from the magnetic susceptibility 

^ ^L5], we just focus on the electronic structure instead of its 



measurements in the HT phase [|14L 
magnetic properties. 

The band structures obtained from LSDA and LSDA+U calculations are shown in Fig. 3. 
Within LSDA (Fig. 3(a)), the bands near the Fermi level (E F ) are mainly derived from V 3d 
states. Since straight V-O-V paths are not present in layered NaV0 2 and instead only nearly 90° 
V-O-V bonds exit, the V 3d states are quite narrow. In the approximately octahedral crystal field, 
the 3d orbitals are split into upper e g and lower t 2g states. As shown in Fig. 3(a), e g derived bands 
range from 1.5 to 2.5 eV and t 2g derived bands lie between -1.5 and 0.5 eV. The splitting between 
t 2g and e g bands is about 1 eV. Under the trigonal crystal field, the t 2g orbitals are further split into 
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one dig and two degenerate e' g orbitals. However, the splitting is much less than the band widths 
so that the t 2g orbitals still cross the E F , which denotes a metallic state within LSDA. That is to 
say, LSDA calculations can not reproduce the insulating nature of NaVC>2 from experiment [14]. 
The LSDA+U scheme|20|] is thus used to count the strong correlation of V 3d electrons. As shown 
in Fig. 3(b), the empty a\ g band is pushed upwards by about 1 eV, and a gap is opened near the 
E F . The system is hence an insulator due to electron correlation and NaV0 2 is indeed a good 
candidate for Mott-Hubbard insulator. 

According to the pure crystal field theory, the a\ g orbital is of lower energy than the e ' g orbitals 
under the trigonal distortion of compression, which is opposite to our LSDA+U results. So it 
is necessary to discuss the controversy on relative order of a\ g -e ' g in such trigonal distortions. 
In Ref. [26], Landron and Lepetit pointed out that this relative order is strongly influenced by 
the e g -e' g hybridization. The e g and e' g orbitals belong to the same irreducible representation 
(E g ) and can thus mix despite the large t 2g -e g energy difference. Such a mix may be small but it 
modulates large energetic factors: the on-site Coulomb repulsions. When the e g -e' g hybridization 
is taken into account, the energy difference AE between the a\ g and e ' g orbitals depends on two 
competitive parts: AE = AE\ + AE 2 = s(a\ g ) - s(e' g ). AE\ includes the kinetic energy, the 
electron-charge interaction, and the interaction with the core electrons. AE 2 denotes the repulsion 
and exchange terms within the 3d shells. Additionally, AE\ and AE 2 both depend on the amplitude 
of the trigonal distortion and are of opposite effect with each other. Under a trigonal distortion 
of compression, if we only consider the crystal field effect (AE{), the a\ g orbital is of lower 
energy than the e' g orbitals (AE < 0). Whereas if we take AE 2 into account, the relative order 
between the a\ a and e' g orbital is reversed (AE > 0), comparing with the crystal field prediction. 
Therefore, LSDA+U calculations predict that the ci\ g orbital is of higher energy than the e ' g orbitals 
in NaV0 2 . In fact, such a controversy has been presented in another compressed triangular system 
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2911 . From the crystal field theory, some authors [27] proposed that the energy 



Na^CoOalaa 

of a\g orbital is lower than the e' g orbitals. However, the LDA+U method[|28[. 



2911 yielded an 



dig orbital of higher energy than the e' g orbitals. Later, the experimental results[30] showed that 
the Fermi surface of the C0O2 layers issues from the d\ g orbital, not at all from the e' g orbitals, 
supporting the LDA+U results. 
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B. IT phase 

The triangular lattice of NaV0 2 exhibits magnetic frustration and spatially anisotropic ex- 
change interactions in the IT phase. As shown in Table II, the G-AFI configuration is the most 
stable state among the five magnetic structures both from LSDA and LSDA+U calculations. By 
a detailed analysis of the magnetic ground state G-AFI, see Fig. 2(a)(i), AFM chains are formed 
along the (110) direction (long V-V bonds), while AFM exchange is also more favorable along the 
(010) direction (short V-V bonds). Considering all the short V-V bonds are completely equivalent, 
both of the (100) and (010) directions should be AFM exchange. Thus, NaV0 2 in the IT phase 
can be regarded as a system with frustration effects. In addition, spatially anisotropic exchange 
interactions may exist in such a triangular spin lattice^], i.e., J\ along the direction of long V-V 
bonds and 7 2 along the two directions of short V-V bonds (Fig. 2(a)(i)). 

In order to describe the magnetic frustration and spatially anisotropic exchange interaction 
more clearly, we estimate the exchange interactions along one of the triangle directions (7i) and 
the other two (7 2 ) (Fig. 2(a)(i)). Since all the configurations exhibit insulator characteristics, the 
spin size of V is stable, and the difference of total energy between C-AFI and G-AFI (C-AFII 
and G-AFII) configurations (See Table II) is so small that the system exhibits a 2D characteristic, 
a nearest neighbor Heisenberg-like Hamiltonian may be a good primary approximation for the 
in-layer magnetic energy. The corresponding 2D spin Hamiltommian can be written as 

H=J 1 ^2s k -S l + J 2 ^2s i -S j (1) 

(k,l) 

where denotes a nearest-neighbor pair (short V-V bond) and (k,l) denotes a next-nearest- 
neighbor pair (long V-V bond). By mapping the obtained total energies for each magnetic state to 
the Heisenberg model, the exchange interactions J% and 7 2 were calculated within this approxima- 
tion: 

2 x (8 x 4 J 2 S 2 ) = E(FM) - E(C - AFII) (2) 

2 x (4 x 4 J 2 S 2 + 4 x 4 J^ 2 ) = E(FM) - E(C - AFI) (3) 

With S=l, we get 7 2 =2.1 meV and 7i=6.1 meV for NaV0 2 in the IT phase, which reflects strong 
spatial anisotropy. The AFM chains are established along the (110) direction (7i) and the inter- 
chain coupling (7 2 ) is frustrated. Moreover, the value of 7 2 /7i=0.3 is so small that this magnetic 



structure can be described as so-called weakly coupled zigzag (S=l) chains model[31]. 
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The integer spins (S=l) are able to weaken the frustration effects in the frustrated systems, as 



in kagome lattice 032Q. Such a lattice has four nearest neighbors with the adjacent triangles on the 



lattice sharing only one lattice point. Interestingly, the triangular lattice can be composed of four 



kagome lattices[27]. Thus, there are some analogous properties in these two frustrated systems. 
Therefore, it is reasonable to suppose that the triangular lattice also has the rule that the half- 
odd-integer spins are more highly frustrated than integer ones. For example, NaTi0 2 (S = |)[|5D 
and LiCr0 2 (S = with half-odd-integer spins are always frustrated even at low T, while 

the magnetic frustration of NaMn0 2 (5=2) is clearly lifted by a structural distortion[2]. Besides 
NaMn02, NaV0 2 is another typical triangular lattice with integer spins (S=l). Therefore, we can 
presume that the magnetic frustration in NaV0 2 can be lifted in some way. 

C. LT phase 

From the discussion above, we expect that NaV0 2 with S=l will show a long-range magnetic 
ordering or a finite ground-state magnetization at low T. As shown in Table II, G-AFII is only 0.3 
meV lower in total energy than C-AFII within LSDA, and both G-AFII and C-AFII have the same 
lowest total energy from LSDA+U results for the LT phase. This reflects the obvious 2D charac- 
teristic of NaV0 2 : the interlayer interaction is much weaker than the intralayer one. As stacking 



antiferromagnetically between layers is observed in the experiment 14], G-AFII state should be 
more favorable at low T. Such a magnetic state denotes the long range 3D magnetic ordering with 
AFM coupled FM chains in V cation layers and interlayer AFM coupling. Obviously, the mag- 
netic frustration is lifted in G-AFII state by the first-order transition at 91.5 K: the lattice distorts 
to another monoclinic (C2/m) with four long and two short V-V bonds reversed, comparing with 
IT phase. 

The lattice distortion, which relieves the frustration, is driven by the formation of 00 in the 
LT phase. Fig. 4 shows the orbital characteristic of V 3d in G-AFII state. Since the orbital 
occupancies of the two inequivalent V (VI and V2 in Fig. 2(a)(ii)) are nearly the same, only the 
density of states (DOS) of VI 3d is shown. The z axis of local coordinate system coincides with 
the V-0 bond of the V0 6 octahedra. In such a coordinate system, d zx and d yz orbitals are mainly 



occupied and d xy orbital is less occupied at al 



V ions. Such an orbital occupancy is consistent 



with the 00 proposed by McQueen and Cava lll4h : the d zx and d yz orbitals are singly occupied 
with all unoccupied d xy orbitals. 
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This 00 relieves the magnetic frustration and stabilizes the long-range magnetic ordering state. 
In view of the weak superexchange interaction resulted from the nearly 90° angle of V-0-V II33ll 
as well as the weak magnetic interaction between adjacent VO2 planes interleaved by a layer of 
Na ions, the V-V direct exchange interaction in-plane should be dominant in such a particular 
crystal structure. Particularly, we only consider the a overlap in V-V direct exchange, which is 
much stronger than the ir overlap. It means that each orbital in a V ion only hybridizes with the 
same orbitals in the two nearest-neighboring V ions. That is to say, d yz orbital hybridizes with 
two neighboring d yz orbitals in the (010) and (010) directions, and d zx orbital hybridizes with two 
neighboring d zx orbitals in the (100) and (100) directions. The repulsions between the occupied 
orbitals (d zx or d yz ) induce the elongation of V-V bonds along four directions: (010), (010), (100) 
and (100). According to the Goodenough-Kanamori (GK) rules [|34fl. the strong AFM coupling 
should exist along these four directions because of the occupation of two orbitals with a overlap. 
At the same time, the less occupancy of d xy orbital leads to V-V bonds contraction as well as a 
weak FM exchange along the (110) and (110) directions. Thus, the four long and two short V-V 
bonds result from bonding via d zx , d yz orbitals, but no bonding of the d xy electrons. In other words, 
such an OO results in the lattice distortion, and consequently relieves the magnetic frustration. 

In the LT phase, another important aspect is that the SOC turns out to be crucial for the small 
magnetic moment of 0.98 [Ib per V 3+ observed experimentally [14]. As shown in Table II, the 
magnetic moments from LSDA+U calculations are much larger than the ones observed in the ex- 
periment. Further to investigate the magnetic moments changing with the particular choice of 
U e jf, we calculate the moments for f/ e //=2-6 eV and find that the moments are not sensitive to 
U e ff\ as shown in Table III, the magnetic moments stay constant within 0.2 ^ as long as the sys- 
tem is an insulator. Since an easy magnetization direction (1 10) is observed in the exprement[14], 
the SOC may play an important role in determining the total magnetic moment. Thus, the SOC is 
included to reinvestigate the magnetic moment. 

Then, we perform LSDA+SOC+U calculations for the favorable magnetic configuration (G- 
AFII), and obtain a local moment of 0.89 fig per V 3+ with 1.65 spin and -0.77 [is orbital 
contributions. This value is half the expected moment (2 //#), but very close to the observed one 
(0.98 hb)- The DOS projected on (2,m) space shown in Fig. 5 reveals the origin of orbital mo- 
ment. Note that the z axis is set to the direction of easy magnetization along (110) now. Since d\ 
and d _i have nearly the same occupancies, the orbital moment only comes from the contribution 
of different occupancies between J 2 and J_ 2 . By further analysis, the J 2 occupancy is less than 
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one half of the <i_ 2 one, which should give an orbital moment between 1 jj, B and 2 ji B . Never- 
theless, there is no surprising that the calculated orbital moment is 0.77 fi B here, because some 
hybridization effects are neglected in above analysis, e.g., the covalence effects with O 2p. Thus, 
the inclusion of SOC leads to a surprising but experimentally sound results. 



IV. CONCLUSIONS 



In Summary, we have investigated the electronic structure and magnetic properties of NaVC"2 
by first-principles calculations. The ti g orbitals are split into upper ci\ g and lower e' g states 
by a trigonal distortion of compression in the HT phase, which is similar to the splitting in 
Na x Co02|28|]. In the IT phase, the crystal symmetry is lowered to C2/m, under which the sys- 
tem behaves like a frustrated spin lattice with spatially anisotropic exchange interactions. Finally, 
a long-range ordered AFM ground state is formed when the magnetic frustration is relieved by 
another lattice distortion resulted from a certain ordering of occupied orbitals at low T. The small 
magnetic moment observed originates from the compensation of orbital moment for the spin mo- 
ment. It is obvious that so many physical phenomena in the triangular lattice are reflected in 
NaV0 2 , suggesting that NaV0 2 is a very good model material for studying 2D triangular lattice 
systems. 



V. ACKNOWLEDGEMENTS 



We thank T. M. McQueen, P. W. Stephens, Q. Huang, T. Klimczuk, F. Ronning, and R. J. Cava 
for providing structural parameters prior to publication. This work was supported by the special 
Funds for Major State Basic Research Project of China(973) under grant No. 2007CB925004, 863 
Project, Knowledge Innovation Program of Chinese Academy of Sciences, and Director Grants of 
CASHIPS, CUHK Direct Grant No. 2060345. Part of the calculations were performed in Center 
for Computational Science of CASHIPS and the Shanghai Supercomputer Center. 



[1] K. Takada, H. Sakurai, E. Takayama-Muromachi, F. Izumi, R. A. Dilanian, and T. Sasali, Nature 
(London) 422, 53 (2003). 



10 



[2] Maud Giot, Laurent C. Chapon, John Androulakis, Mark A. Green, Paolo G. Radaelli, and Alexandras 

Lappas, Phys. Rev. Lett. 99, 247211 (2007). 
[3] C. Darie, P. Bordet, S. de Brion, M. Holzapfel, O. Isnard, A. Lecchi, J. E. Lorenzo, and E. Suard, Eur. 

Phys. J. B 43, 159 (2005). 

[4] F. Reynaud, D. Mertz, F. Celestini, J.-M. Debierre, A. M. Ghorayeb, P. Simon, A. Stepanov, J. Voiron, 

and C. Delmas, Phys. Rev. Lett. 86, 3638 (2001). 
[5] S. Yu. Ezhov, V. I. Anisimov, H. F. Pen, D. I. Khomskii and G. A. Sawatzky, Europhys. Lett. 44, 491 

(1998). 

[6] 1. 1. Mazin, Phys. Rev. B 75, 094407 (2007). 

[7] Albert J. W. Reitsma, Louis Felix Feiner and Andrzej M. Oles, New J. Phys. 7, 121 (2005). 
[8] F. Mila, F. Vernay, A. Ralko, F. Becca, P. Fazekas and K. Penc, J. Phys.: Condens. Matter 19, 145201 
(2007). 

[9] J.-H. Chung, Th. Proffen, S. Shamoto, A. M. Ghorayeb, L. Croguennec, W. Tian, B. C. Sales, R. Jin, 
D. Mandrus and T. Egami, Phys. Rev. B 71, 064410 (2005). 
[10] Y. Tokura and N. Nagaosa, science, 288, 462 (2000) 

[1 1] J.-Q. Yan, J.-S. Zhou, and J. B. Goodenough, Phys. Rev. Lett. 93, 235901 (2004). 
[12] Peter Horsch, Andrzej M. Oles, Louis Felix Feiner, and Giniyat Khaliullin, Phys. Rev. Lett. 100, 
167205(2008). 

[13] H. F. Pen, J. van den Brink, D. I. Khomskii, and G. A. Sawatzky , Phys. Rev. Lett. 78, 1323 (1997). 
[14] T. M. McQueen, P. W. Stephens, Q. Huang, T. Klimczuk, F. Ronning, and R. J. Cava, Phys. Rev. Lett. 

101, 166402 (2008). 
[15] M. Onoda, J. Phys. Condens. Matter 20, 145205 (2008). 

[16] H. F. Pen, L. H. Tjeng, E. Pellegrin, F. M. F. de Groot, G. A. Sawatzky, M. A. van Veenendaal and C. 

T. Chen, Phys. Rev. B 55, 15500 (1997). 
[17] The difference between 93 K and 91.5 K has to do with the sample being measured on warming 

or on cooling. The IT to LT transition temperature is T=91.5 K here, corresponding to the lattice 



parameters measured on cooling. And the slight lattice-parameter differences from Ref. Ill4ll are within 

experimental error. 
[18] P. Blaha et al., http: // www. wien2k. at. 
[19] John P. Perdew and Yue Wang, Phys. Rev. B 45, 13244 (1992). 

[20] Vladimir I. Anisimov, Jan Zaanen, and Ole K. Andersen, Phys. Rev. B 44, 943 (1991). 



11 



[21] S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys and A. P. Sutton, Phys. Rev. B 57, 1505 
(1998). 

[22] G. K. H. Madsen and P. Novak, Europhys. Lett. 69, 777 (2005). 

[23] V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. Czyiyk and G. A. Sawatzky, Phys. Rev. B 48, 
16929 (1993). 

[24] Robert Laskowski, Peter Blaha, and Karlheinz Schwarz, Phys. Rev. B 67, 075102 (2003). 

[25] We performed LSDA+U calculations for the two most stable configurations (G-AFI and G-AFII) both 

in the IT and LT phases with U eff =2, 3, 5, 6 eV. The results show that the G-AFI (G-AFII) is always 

the most stable state in the IT (LT) phase in such U e f / range. 
[26] Sylvain Landron and Marie-Bernadette Lepetit, Phys. Rev. B 77, 125106 (2008). 
[27] W. Koshibae and S. Maekawa, Phys. Rev. Lett. 91, 257003 (2003). 
[28] Sylvain Landron and Marie-Bernadette Lepetit, Phys. Rev. B 74, 184507 (2006). 
[29] Liang-Jian Zou, J.-L. Wang, and Z. Zeng, Phys. Rev. B 69, 132505 (2004). 

[30] M. Z. Hasan, Y.-D. Chuang, D. Qian, Y. W. Li, Y. Kong, A. P. Kuprin, A.V. Fedorov, R. Kimmerling, 
E. Rotenberg, K. Rossnagel, Z. Hussain, H. Koh, N. S. Rogado, M. L. Foo, and R. J. Cava, Phys. Rev. 
Lett. 92, 246402 (2004). 

[31] Zheng Weihong, Ross H. McKenzie and Rajiv R. P. Singh, Phys. Rev. B 59, 14367 (1999). 

[32] Swapan K. Pati and C. N. R. Rao, J. Chem. Phys. 123, 234703 (2005). 

[33] J. B. Goodenough, Phys. Rev. 117, 1442 (1960). 

[34] J. B. Goodenough, Magnetism and the Chemical Bond (Interscience Publishers, New York, 1963); J. 
Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 



12 



TABLE I: The lattice parameters of NaV0 2 at 100, 91.5 and 20 K 



Temperature 


100 K 


91.5 K 


20 K 


space group 


R-3m 


C2/m 


C2/m 


a(A) 


2.9959(1) 


5.1758(1) 


5.2223(2) 


b(A) 


2.9959(1) 


3.0078(1) 


2.9755(1) 


c(A) 


16.0996(1) 


5.6340(1) 


5.6492(3) 


a(°) 


90 


90 


90 


PC) 


90 


107.629(1) 


108.335(1) 


7(°) 


120 


90 


90 


x,y,z(Na) 


3a(0, 0, 0) 


2a(0, 0, 0) 


2a(0, 0, 0) 


x,y,z(V) 


3b(0, 0, 0.5) 


2d(0, 0.5, 0.5) 


2d(0, 0.5, 0.5) 


x,y,z(0) 


6c(0, 0, 0.2339(0)) 


4i(0.2368(7), 0, 0.6989(5)) 


4i(0.2296(5), 0, 0.7005(4)) 
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TABLE II: The total energy E (meV/ (8f. u.)), magnetic moment M (/zb) per V 3+ and band gap E g (eV) in 
different magnetic states. 





Configuration 


FM 


C-AFI 


G-AFI 


C-AFII 


G-AFII 


IT (T=91.5 K) 


E 


514 


4 





128 


132 


LSDA 


M 


1.63 


±1.37 


±1.37 


±1.34 


±1.34 


LT (T=20 K) 


E 


574 


152 


156 


4 







M 


1.51 


±1.35 


±1.35 


±1.38 


±1.38 




E 


329 


65 





192 


189 


IT (T=91.5 K) 


M 


1.71 


±1.65 


±1.65 


±1.66 


±1.66 


LSDA+U 


Eg 


1.1 


1.2 


1.1 


1.4 


1.5 




E 


395 


286 


288 








LT (T=20 K) 


M 


1.71 


±1.66 


±1.66 


±1.65 


±1.65 




Eg 


1.2 


1.4 


1.4 


1.1 


1.1 
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TABLE III: The band gap E g (eV) and magnetic moment M{hb) per V 3+ for different U e f f (eV) in G-AFII 



configuration. 



Ueff 


2 


3 


3.6 


5 


6 


E 9 


0.02 


0.9 


1.1 


1.8 


2.3 


M 


±1.57 


±1.63 


±1.65 


±1.69 


±1.71 
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FIG. 1: (color online), (a) The compressed octahedron of VO2 layers. The z axis is the three fold axis of the 
VC>6 octahedron, a represents the O-V-0 angle, (b) a angles in the HT, IT and LT phases. 

FIG. 2: (a) Two different magnetic patterns in V-V plane (i) I-type AF and (ii) II-type AF. The solid (dashed) 
lines represent short V-V bonds in the IT (LT) phase and dashed (solid) lines represent long V-V bonds in 
the IT (LT) phase. J\ (J2) denotes exchange interaction along the direction of dashed (solid) lines, (b) 
Schematic representation of five magnetic configurations used in our calculation. Only V atoms are drawn. 
Filled (open) circles indicate spin up (down) moments. 

FIG. 3: The spin-majority band structure of NaV0 2 from (a) LSDA and (b) LSDA+U (U eff =3.6 eV) 
calculations for the fixed structure at 100 K. The band with (a) 3d and (b) a\ g character is marked. 

FIG. 4: Density of states (DOS) of NaV0 2 calculated by LSDA+U (U eff =3.6 eV) in the G-AFII state for 
the fixed structure at 20 K. Besides total 3d state, all the d zx , d yz and d xy orbitals in local coordinate system 
for VI (Fig. 2(a)(i)) ion are depicted. Solid (dashed) lines denote the spin-up (down) states. 

FIG. 5: Density of states (DOS) of NaV0 2 calculated by LSDA+U+SOC (t/ e //=3.6 eV) in the G-AFII state 
for the fixed structure at 20 K. Solid (dashed) lines denote the spin-up (down) states. 
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